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Abstract 

We study how the environment of photonic band gap (PBG) materials affects entanglement 
dynamics of qubits. Entanglement between the single qubit and the PBG environment is investi- 
gated through the von Neumann entropy while that for two initially entangled qubits in this PBG 
reservoir is through concurrence. Dynamics of these measurements are solved in use of the frac- 
tional calculus which has been shown appropriate for the systems with non-Markovian dynamics. 
Entropy dynamics of the single qubit system reveals that the coupling with the PBG reservoir 
prevents decoherence of the qubit through the steady entropy with non-zero value. The effect of 
PBG reservoir on the concurrence of the two-qubit system leads to the long-time entanglement 
preservation. The concurrence dynamics shows that unphysical entanglement trapping does not 
exist in the system with the qubit frequency lying outside the PBG region. Long-time memory 
effect of the PBG reservoir occurs only for the qubit frequency in the PBG region. Entanglement 
mechanisms resulting from this long-time memory effect are discussed. 
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I. INTRODUCTION 



Photonic band gap (PBG) materials, in which no travehng electromagnetic (EM) modes 
are allowed in the PBG of their band structures, enrich the content of quantum optics and 
quantum information processing through its non-Markovian dynamics. Periodic structures 
constructing these PBG materials posses the controllable features by the lattice parameters. 
The development of these PBG materials has attracted the interest over the past twenty 
years starting from Yablonovitch's seminar paper |Tj, where the periodic structures with 
dielectric variation of photonic crystals (PhCs) were proposed to be a means to inhibit 
spontaneous emission (SE). Shortly thereafter, John el al. found that the PhC mediates 
a coherent interaction between the emitter and the previously emitted radiation [2]. The 
memory effect in this coherent interaction leads to non-Markovian emission dynamics such as 
non-exponential decay and Rabi-like oscillation in the spectral vicinity of a photonic band 
edge |3]. This non-Markovian emission dynamics enriches the content of quantum optics 
and quantum information processing in the PBG materials. Recently, tunable artificial 
crystals consisting of Josephson junction arrays have been designed as a PBG structure 
with tunable frequency range through the applied external flux [1H6]. This PBG structure 
of superconducting circuits with more adjustable parameters than PhCs might ultimately 
form the primitive blocks of quantum computers. 

When two initially entangled qubits are placed in the PBG materials, it has been reported 
that the controllable SE of the individual qubit can be used to manipulate the entanglement 
of the two qubits. For example, Bellomo et al. [7j suggested high value of entanglement trap- 
ping between two qubits can be achieved by embedding in an isotropic PBG environment. 
The direct link between the entanglement dynamics and the controllable SE dynamics of the 
single qubit is the basis of achieving entanglement preservation. Besides, Amri et al. [8] in 
exploring the entanglement dynamics of two qubits embedded in isotropic and anisotropic 
PhCs found that the initially entangled qubits will become disentangled as the excited qubits 
spontaneously emit the photons into the environment. Whenever the individual qubit forms 
the bound state with the emitted photon, the initial entanglement between the two qubits 
can be partially recovered. Deep understanding about what the entanglement mechanism 
is between the qubits by the PBG reservoir through the non-Markovian dynamics of the 
qubits is now indispensable due to the possibilities for realizing strong-coupling conditions 
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in the exploration systems of quantum information processing [9l [TO] . 

Here we consider the system of qubits placed in a PEG material such as a PhC or tunable 
artificial crystal with directional-dependent dispersion relation near the photonic band edge. 
This anisotropic band structure of the PEG reservoir can be described by the effective- mass 
approximation [HI [12] and leads to the photon density of state (DOS) of the reservoir 
proportional to ^uj — ujc^iuj — Wc), where ujc is the band edge frequency and 0(w — Wc) 
is the Heaviside step function. Since the traveling EM modes are forbidden in the PEG 
region, the emitted photon from the qubit can't travel through the PEG material and thus 
is localized around the qubit to form a bound state. This bound state carries the memory of 
the PEG reservoir on the excited-state qubit. This memory effect controls the dynamics of 
the qubit through the delay Green function or memory kernel. The fractional calculus [13] 
has been shown an appropriate mathematical method for solving these physical systems with 
memory effect. Using the fractional calculus to investigate the non-Markovian SE dynamics 
of the PhCs we obtained the results consistent with the experimental observation [121 [H]. 

Here we use the fractional calculus to investigate how the entanglement dynamics of the 
qubits are affected by the memory effect of the PEG reservoir. The systems we consider 
include one single qubit and two entangled qubits surrounded by the PEG materials, respec- 
tively. The entanglement between the single qubit and the PEG reservoir is investigated 
through the qubit 's excited-state probability and von Neumann entropy while the entan- 
glement between the two qubits is through the concurrence. The probability dynamics of 
the single qubit exhibits non-Markovian and Markovian behavior as the qubit frequency lies 
inside and outside the PEG region, respectively. The memory effect of the PEG reservoir 
on the excited-state qubit leads to the non-Markovian dynamics of the qubit through the 
formation of photon bound state. This memory effect leads to the entanglement preservation 
between the single qubit and the PEG reservoir through the steady entropy with non-zero 
value. In the two-qubit system, we find that unphysical entanglement trapping does not 
exist in the system with qubit frequency lying outside the PEG region. These new results of 
the PEG on the entanglement dynamics are different from those of previous studies which 
predicted the existence of entanglement trapping in the system with qubit frequency outside 
the PEG region [TJ |8]. The accuracy of these results is certified by the use of appropriate 
mathematical method of fractional calculus for non-Markovian dynamics and the reason- 
able results of physical phenomenon. Eesides, we discuss the entanglement mechanisms for 
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the two-qubit system with different prepared states. PBG reservoirs provide the qubits a 
better environment for preserving entanglement through interchanging the trapped photon. 
These accurate results about the effect of PBG reservoir on the qubit's dynamics can be ap- 
plied directly to the experimental systems of tunable artificial crystals with PBG structure 
through the circuit quantum electrodynamics (circuit QED) setup. These artificial crystals, 
composed of periodic arrays of circuit elements of Josephson junction, have been designed 
to have a gap in the EM spectrum through the array parameters [lH6l [T514T8] . As qubits 
are placed in the middle of these arrays of artificial crystals, the emitted photons from the 
qubits in the PBG environment is equivalent to that of qubits embedded inside a PhC. 

The paper is organized as follows. In Sec. II, starting from the quantum theory of one 
single qubit in an anisotropic PhC, we express the kinetic equation of this interaction system 
as a fractional Langevin equation in use of the fractional calculus. By solving the fractional 
Langevin equation analytically, we express the time evolution of the qubit's quantum state 
as a reduced density matrix. Dynamics of probability and von Neumann entropy associated 
with the matrix elements are discussed. In Sec. Ill, we construct the reduced density 
matrix of the two-qubit system through the single-qubit matrix and the two-qubit basis. The 
entanglement dynamics of the two-qubit system is discussed on the basis of the concurrences 
of two initially entangled states. The entanglement mechanisms for the non-Markovian 
systems with different prepared states are revealed in this section. Finally, we summarize 
our results in Sec. IV. 

II. PROBABILITY AND VON NEUMANN ENTROPY OF ONE SINGLE QUBIT 

As a single qubit is surrounded by an anisotropic PBG material with cut-off photonic 
DOS shown in Fig. 1, the total Hamiltonian of this interaction system can be expressed as 

H = lhw,oi\l) (1| - |0) (0|) + X] K^ata,^ + thJ2 9^ (1| ® 4 + h.c), (1) 

k k 

where |1) and |0) stand for the excited and ground states of the qubit with transition 
frequency uio nearly resonant with the frequency range of PBG; and at are the anni- 
hilation and creation operators of photon with wavevector k and corresponding frequency 
in the reservoir; coupling strength between the qubit and the k photon is specified by 
9k ~ [ 2£oa)-y ] ^^fc ' ' '^hich is assumed to be independent of the qubit's position with 
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fixed qubit's dipole moment dio — dioUd- Here the symbols V, eg and eo stand for the sample 

— * 

volume, polarization unit vector of the k photon and the dielectric constant, respectively. 

If we use the coordinate {6, 0) on the Bloch sphere to parameterize the state of the qubit, 
the initial state can be written as 



1^(0)) 



(2) 



with the initial coordinate (^o; 0o) = (O, 0) for the excited state and (tt, 0) for the ground 
state. Here the photon vacuum state and one- photon state are expressed as |0j^) and 
respectively. We assume Cj^(O) = for no initial correlation between the qubit and PEG 
reservoir. This system evolves with time through 



u,{t)e'^''cos{^-^)\l) +ua{t)sin{^-^)\0) 



with initial condition Up{0) — 1, Ud{0) — 1 and Cg(0) = 0. Here Up{t) and Ud{t) stand for the 
excited-state and ground-state probability amplitudes of the qubit with the photon vacuum 
state while Cg(t) for the qubit in its ground state with one photon in the reservoir. 

Substituting Eq. (3) into the time-dependent Schrodinger equation then taking inner 
product with {ip{t)\, we obtain the coupled amplitude equations as 

d 



-uJt) ~ 

dt ' cos 



d 
dt 



Ckit) = g^up{t)cos{0o/2)e 



iClrt 



and 



d 
dt 



Ud{t) = 



(4) 
(5) 
(6) 



with the detuning frequency Qv 



The last equation of motion yields Ud{t) = Ud{0) = 1 meaning that the ground-state 
probability amplitude will not evolve with time. Integrating Eq. (5) and substituting into 
Eq. (4), we have 



JtMt)--l G{t - T)Up{T)dT (7) 

with the memory kernel G{t — r) = X^j^ 5'|e~*'^**-*~^'*. This equation reveals that the future 
state of the qubit is related to the memory of the reservoir in its previous state through the 
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memory kernel. As the qubit is put in free space, the memory kernel has the form of a Dirac 
delta function G{t — r) oc 6{t — r) corresponding to continuous photon DOS. In this case, 
the reservoir responses only at an instant time r = t which leads to the excited amplitude of 
the qubit decaying exponentially with time. This Markovian result in free space manifests 
that the qubit loses all memory of its past and decays quickly to its ground state. 

For the anisotropic PBG reservoir discussed here, the memory kernel manifests its mem- 
ory effect within the entire interval (0,t) through G{t — r) = ^JJ^yii^ e~''[^^/^~'^(*~'^)] with 
the coupling constant /S"*^/^ = (co'^Q(i^Q^/u^)/(167reo^c^) and the detuning b = Uio — Uc of 
the qubit frequency Uiq from the band edge frequency Uc [12] ■ This memory kernel is de- 
rived from the anisotropic dispersion relation of the PBG material. Near the band edge 
frequency Uc, the anisotropic dispersion relation has a vector form and can be expressed 
by the effective-mass approximation as [H] uj: ~ Uc + A (jc — with the curvature 
A = fuc/kl = fc^/oJc signifying its directional-dependent values through the scaling factor 
/. This dispersion relation leads to the memory kernel expressed by the cut-off photon DOS 
= 4^\/^©(^-^c) asG(t-r) = '^Mji duj^e''^'^-^''^^^'-^\ Substituting this 
memory kernel for the anisotropic PBG reservoir into Eq. (7), we obtain 

d , , /3V2ei3-/4 rt i^ym-r) 



This equation can be further simplified by making the transformation Up{t) = e^^^Upit) to 
give 

jUM^mM-^-^lj^^dr. (9) 



Conventionally, this integro-differential equation is solved through Laplace transform which 
leads to the fractal phenomenon of the system for the memory kernel in Laplace image |T9] . 
This fractal phenomenon would result in the stochastic nature of the dynamical behavior of 
the system which appears as the non-Markovian dynamics. This non-Markovian dynamics 
can be solved analytically by the fractional calculus. 

Here we solve the non-Markovian dynamics of the single qubit system by the fractional 
calculus. The detailed derivation can be found in Appendix A. We obtain the analytical 
solution of Eq. (9) as 

U,{t) = X [n'^*(l/2, Y^) - YiE,{l/2, Yi) + Y.e"^?' - y^e^^^*] , (10) 
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for PIP ^ 5; and 



Up{t) = -2 



^3/2gi37r/4 



J3/2 



(11) 



ior p/p — S. The time evolution of the system is thus analytically expressed as 



h) + Y.^kim^\W (12) 



with Up{t) = e^'%{t), ua{t) = uM = 1 and Y.k = 1 - [u^{t)cos{^)f - [s^n(|)]^ 

The reduced density matrix of a single qubit in the anisotropic PBG reservoir can also be 
expressed analytically as 



\up{t)\^ cos^{^) lu;{t)e-'^°sin{eo) 
\up{t)e^'^°sin{9o) 1 - 1^^(^)1^0052(1) 



Pii(^) Pio(^) 
Poi(^) Poo(^) 



(13) 



le-^'t'osin{9o) sin^{^) j 
We show the effect of PBG reservoir on the excited-state probability and von Neumann 

entropy of the single qubit in Fig. 2. For the qubit frequency uio hes above the edge 
frequency ouc or the positive detuning case [S/P — {ujiq — uJc)/P > 0], the excited-state 
probability exhibits Markovian behavior with exponential decay. The initially excited qubit 
loses almost its energy into the PBG reservoir and decays to the ground state. During this 
decaying process, the qubit loses all memory of its past through emitting a photon. For 
the negative detuning case with the qubit frequency lying inside the PBG region [S/f3 — 
{ouio — ijJc)/P < 0], the system manifests decay and oscillatory behavior before reaching 
a non-zero steady-state value. This non-Markovian dynamics reveals that the qubit loses 
partial of its energy in the very beginning period of time and then preserves the remaining 
energy in the end. In this case, the emitted photon from the qubit 'sees' no allowed DOS 
in the PBG region so that it is localized around the qubit through the reflecting and re- 
absorbed processes and forms the bound state with the qubit. The formation of the bound 
state preserves the remaining energy of the qubit and the memory of the PBG reservoir on 
the excited qubit. This memory effect leads to the distinct dynamical behavior of the non- 
Markovian system from the Markovian one. When a quantum operation with multiple stages 
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with the initial one p^(0) 



is performed on the qubit, the system with memory effect will experience non-Mar kovian 
processes of operation with the pre-stage operation consecutively affecting the following one. 
On the other hand, the qubit in the Markovian system without memory effect will experience 
stochastic processes of operation with each stage operation independent of the other. 

In Fig. 2(b), we show the dynamics of von Neumann entropy in terms of the eigenvalues 
of the reduced density matrix in Eq. (13). This quantum entropy is defined as S(t) = 
—Tr [p'^{t)logp^{t)] = — J2i where Aj are the eigenvalues of the density matrix p*(t). 

The entropy measures the correlation between the reservoir and the qubit through quantify- 
ing the amount of the information we will gain after we measure a quantum state. The larger 
the entropy value is, the stronger correlation exists. For the density matrix in Eq. (13), the 
eigenvalues can be obtained easily as A± = ^ |l ± y^'l — 4cos^ (y) [|'Up(t)|^ — |'Up(t)|^] j>. 
We plot the entropy dynamics in Fig. 2(b) based on these eigenvalues. This quantum en- 
tropy evolves with time because the correlation with the PBG reservoir leads to the amount 
of information varying with time. For the initial pure state of the excited qubit, the entropy 
exhibits its minimal value of zero at t = and reaches its maximal value log2 = 0.693 at the 
very beginning of time. After a period of time of the decay timescale, the entropy becomes 
steady with non-zero value in the non-Mar kovian system {6/(3 < 0) and zero in the Marko- 
vian system (5//3 = 2). These results show that the initial pure system becomes maximally 
mixed in the very beginning period of time because of correlation with the reservoir. As 
the qubit equilibrates with the PBG reservoir, the system becomes less mixed. Contrarily, 
in the Markovian system, the qubit decays to its another pure state (ground state), as it 
quickly equilibrates with the reservoir. In this case, the stored information in the qubit 
is lost because of the correlation with the Markovian reservoir [20j. However, in the non- 
Markovian system, the strong correlation with the PBG reservoir leads to the preservation 
of the mixed state. The larger non-zero steady entropy value of the non-Markovian system 
has, the greater amount of information we will gain upon measuring this mixed state. Since 
the system with qubit frequency deeper inside PBG region has the larger non-zero steady 
value of entropy, we can use the absolute value of qubit's detuning frequency \6\ / (3 to char- 
acterize the degree of the non-Markovian effect from PBG reservoir. Correlation with the 
PBG reservoir prevents losing stored information from the single qubit through the large 
value of the steady entropy. 
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III. ENTANGLEMENT OF TWO QUBITS 

In the following we consider a system consisting of two independent parts, each with 
a qubit coupled to the PBG reservoir with anisotropic one-band model. The two qubits 
A and B with the same transition frequency Uio are assumed to be identical and initially 
entangled. The distance between the two qubits is larger than the spatial correlation length 
between the individual qubit and PBG reservoir. For each part of the subsystem with one 
qubit and the corresponding PBG reservoir, the reduced density matrix for the single qubit 
can be described by Eq. (13) with the superscript s = A,B. For the two-qubit system, the 
reduced density matrix p{t) can be constructed through the single-qubit matrix p'^(t) and 
the two-qubit basis |11) = |1) , |10) = |2) , |01) = |3) , |00) = |4) [SET]. This reduced density 
matrix for the two-qubit system is a Hermitian matrix with its elements Pij{t), i,j = l 4, 
depending only on the initially entangled condition and the excited-state amplitude Up{t) of 
the single qubit. For the initially entangled states |$) = a |01) + 7 |10) = a |3) + 7 |2) and 
= a |00)-|-7 |11) = a |4)-|-7 |1) with a being real, 7 = I7I e**^ and q;^-|-|7|^ = 1, the matrix 

4 2 r 2 ~l 2 

elements are pii(t) = pii(O) \up{t)\ , p22(i) = Pii(O) \up{t)\ [1 - \up{t)\ \ + p22(0) \up{t)\ , 

P33(t) = Pll(O) \Up{t)f [1 - \Up{t)f] +P33(0) \Up{t)\\ piiit) = 1 - [pil(t) +P22(t) +P33(t)], 

Piiit) = pi4(0) \up{t)f and p23{t) = P23(0) \up{t)f. Here the initial elements are Pij(O) = 
(1$) and (1^) {"^Dij for the prepared |$) and states. 

In order to investigate how the PBG reservoir affects the entanglement of the two qubits, 
we calculate the concurrence of the system which provides the direct measurement of the 
entanglement [22]. For the initially entangled states |$) and the concurrences are given 
by [71E3] 

C$(t) = max. jo, 2 |p23(t)| - 2^pu{t)pM{t)] (14) 

and 

C,t(t) = max. 1 0,2 |pi4(t)| -2 v/p22(t)p33(t)}- (15) 
Substituting the above matrix elements into the concurrences, we arrive at 

C^{t) = max. |o,2aVl - \up{t)f^ (16) 

and 

C^{t) = max. jo, 2VI - \up{t)f |a - Vl - [l - |Mp(t)|^]}}. (17) 
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These expressions of the concurrences have the neat forms of C$(t) = \up(t)f and C^(t) = 
|Mp(t)|'^ if the two qubits are maximally entangled through the Bell states |$)^ = (|01) ± 
\10))/V2 and 1^)^ = (|00) ± |ll))/v^. 

In Fig. 3, we show the concurrence dynamics of the two qubits varying with the initial 
entangled degree of the prepared states |$) and |\E'). For the prepared |$) state, the 
concurrences with qubit frequency lying inside and outside PBG are respectively displayed 
in Fig. 3(a) (non-Markovian system, S/P = —5) and Fig. 3(b) (Markovian system, 6/(3 = 2) 
while those for the prepared state are in Fig. 3(c) (non-Markovian system, = —5) 
and Fig. 3(d) (Markovian system, 6/(3 = 2). The non-Markovian system for the prepared |$) 
state in Fig. 3(a) exhibits promising entanglement trapping for all initial entangled degree 
a ^ 0. This entanglement trapping phenomenon is linked to the excited-state probability 
P{t) = \up{t)f of the single qubit in steady bound state with non-zero value. This behavior 
induced by the memory effect of the PBG reservoir on the excited single qubit is evidently 
different from that of the Markovian system in Fig. 3(b), where C$(t) decays exponentially 
and vanishes permanently. The backaction effect of the PBG reservoir on the non-Markovian 
system leads to the memory effect of the single qubit system and the entanglement trapping 
phenomenon in the two-qubit system. Besides the backaction effect of the PBG reservoir, 
the dependence of the entanglement dynamics on the initially entangled degree is also 
clearly visible in Fig. 3(a). The stronger the initial entanglement is, the larger steady 
value of the concurrence reaches. It thus can be deduced from Fig. 3(a) that the maximally 
entangled |$) state occurs at = 1/2 which leads to the Bell states |$)^ = (|10)±|01))/v^. 
Entanglement prefers to exist in the system with the two qubits sharing just one excitation 
|24j . This fact of Bell state being maximally entangled |$) state is also revealed in the 
Markovian system of Fig. 3(b). The initial entanglement of this Markovian system is 
followed by a period of complete disentanglement. The period of time it takes to reach 
the disentanglement depends on the degree of the initial entanglement. This time period is 
the longest as the initial entangled degree equal to 1/2. The fact that the system with 
two qubits sharing just one excitation exhibits maximally entangled reveals that the two 
qubits preserve their entanglement through interchanging the trapped photon in the PBG 
reservoir. Since the strength of this interchange is associated with the relative position of 
the two qubits, we believe the entanglement in the PBG system depends strongly on the 
position of the two qubits. This will be our future topic to discuss how the relative position 
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of the qubits in the PBG material affects the entanglement of the system. 

Entanglement dynamics for the prepared |\E') state in Fig. 3(c) exhibits the similar 
behavior to those in Fig. 3(a) where entanglement trapping occurs for o? > 1/3. However, 
for < 1/3, entanglement of the non-Markovian system vanishes permanently after a short 
period of time, similar to the Markovian case in Fig. 3(d) j25l [26]. For the prepared |\E') 
state in Figs. 3(c) and (d), the maximally entangled state occurs at = 0.88. The system 
is maximally entangled when the two qubits are almost in the same state coupled with two 
trapped photons. This fact reveals that the system preserves the entanglement through the 
dipole-dipole interaction (DDI) between the two qubits. The DDI has large strength when 
the dipoles of the two qubits have the same orientations. The stronger DDI between the 
two qubits, the better entanglement trapping the system exhibits. 

When we further compare the non-Markovian systems in Fig. 3(a) and (c), we find that 
the dynamical difference between these two systems stems from the different entanglement 
mechanisms. PBG reservoir provides the qubits a better environment for preserving entan- 
glement through interchanging the trapped photon than through DDI. This fact is more 
evident in Fig. 4(a) and (b) where we show the concurrence dynamics for the prepared Bell 
states 1$)^ = (|10)±|10))/v/2 and |^)^ = (|11) ± |00))/y2. The steady values of the entan- 
glement trapping for the prepared Bell state |$) are larger than those for the |\E') state with 
the same detuning frequency 6/p. This fact implies that the product of the probability of the 
states with only one excited qubit (|10), |01)) is larger than that of the totally excited state 
(1 11)) and ground state (|00)) of the two-qubit system. The larger product of probability of 
the states with only one excited qubit corresponds to the greater strength of interchanging 
the trapped photon than strength of DDI. The same result can also be observed in Fig. 
4(c) where we show the steady values of the concurrences for the two prepared Bell states. 
For the same qubit's detuning frequency 6/(3, the system with prepared |$)^ state exhibits 
larger steady values than those of the I^E')^ state which results from the greater strength of 
interchanging the trapped photon. 

When we consider the different detuning frequencies of the qubits in Fig. 4(c), the system 
displays the entanglement trapping only for the frequency lying within the PBG. That is, 
there does not exist entanglement trapping in the system with the qubit frequency lying 
outside PBG region. The strength of entanglement trapping can be indexed by the absolute 
value of detuning frequency \S\ / (3 inside the PBG region. The deeper the qubit frequency 
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is detuned inside the PBG region, the stronger entanglement trapping yields. 

IV. CONCLUSION 

We have investigated the effect of PBG on the dynamics of qubits surrounded by the PBG 
reservoir. In use of the fractional calculus, we solve the qubit's dynamics and analytically 
express the excited-state probability and von Neumann entropy for a single qubit and the 
concurrence for the two-qubit system. The system with qubit frequency lying inside PBG 
exhibits strong non-Markovian nature due to the formation of photon bound state with 
the qubit. This strong non-Markovian nature leads to the long-time memory effect of the 
PBG reservoir on the excited-state qubit. Systems with qubit frequency lying inside the 
PBG experience non-Markovian processes of operation while those outside the PBG will 
experience stochastic processes. The deeper inside the PBG the system is detuned, the 
higher degree of non-Markovianity the system exhibits. The long-time memory effect leading 
to this non-Markovianity leads to the preservation of entanglement between the single qubit 
and the PBG reservoir through the steady entropy with non-zero value. The large value of 
the steady entropy ensures that the quantum information can be stored in the qubit for a 
long period of time. This long-time memory effect disappears as the qubit frequency lying 
outside the PBG region. 

As the two-qubit system is concerned, the results of concurrence dynamics reveal how the 
effect of PBG reservoir controls the entanglement through the initial prepared conditions 
and the degree of non-Markovianity. The entanglement mechanism for the prepared states 
1$) is interchanging the trapped photon and while for the prepared |^) state is dipole-dipole 
interaction. The PBG reservoir provides the qubits a better environment for preserving en- 
tanglement through interchanging the trapped photon than through dipole-dipole interac- 
tion. How the degree of non-Markovianity affects the entanglement is displayed in the steady 
values of concurrence with respect to the detuning frequency of the qubit. The deeper the 
qubit frequency is detuned into the PBG, the higher degree of the non-Markovianity and the 
stronger entanglement trapping yields. The system exhibits the non-Markovian nature of 
achieving entanglement only for the qubit frequency located inside the PBG. Unphysical en- 
tanglement trapping does not exist in the system with the qubit frequency lying outside the 
PBG region. These new predictions of the entropy and concurrence dynamics are different 
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from those of the previous theoretical studies which predicted the existence of entanglement 
trapping in the system with qubit frequency lying outside the PBG region. The accuracy 
of these predictions is based on the appropriate mathematical method of fractional calculus 
for non-Markovian dynamics and reasonable results of physical phenomenon. The results 
about how the PBG reservoir affects the qubit's dynamics can be applied directly to the 
experimental systems of artificial crystals with PBG structured reservoir. We believe the 
entanglement mechanisms between the qubits and PBG reservoir hold the key technology 
of applying the PBG materials to the primitive blocks of quantum computers. 
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Appendix A: CALCULATION OF NON-MARKOVIAN DYNAMICS OF SINGLE 
QUBIT SYSTEM THROUGH FRACTIONAL CALCULUS 

The equation of motion of this qubit system in simphfied form is given in Eq. (9). We 
use the fractional time derivative, one of the operators of fractional calculus, to express the 
integral term of this kinetic equation. In the well-known Riemann-Liouvile definition, the 
fractional time derivative operator d^/df^ is expressed as 



for n < u < n + 1 and T{x) being the Gamma function. The integral term of Eq. (9) 
/o (f-r^)3/2 '^'T can be expressed as this fractional operator with order u = 1/2 and n = 
0. Through this expression and applying the integral operator (d^^/dt^^) and fractional 




(Al) 
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differentiation operator d}/'^ /di}/'^ to Eq. (9), the fractional kinetic equation is obtained as 

jl/2 J-1/2 2/?V2eW4 f-1/2 

+ «^C,W + = (A2) 

This kinetic equation governing the future of the qubit indicates a subordinated stochastic 
process directing to a stable probability distribution because of its interacting with the 
PBG reservoir. This fractional differential equation can be solved analytically through the 
fractional Laplace transform. The basic formula used here is 

with the Laplace variable s. By performing the fractional Laplace transform on Eq. (9), we 
can obtain 



^P^^) s + i(5 + 2^V2eW4si/2 ip/2 
where Up{s) is the Laplace transform of Up{t). As we express the denominator part of this 
Laplace transform as a sum of partial fractions through the roots of the indicial equation 
y^+2/3-'^/^e*'^/'^y//^/^+i5 = with the variable Y converted from s^/^, the analytical solution 
of the Eq. (A2) can be obtained by means of the inverse Laplace transform. There exist two 
kinds of roots in the indicial equation: degenerate root and non-degenerate ones. For the 
degenerate case with — 5, the partial-fractional form of Up{s) — 1/ (^y/s -\- ^-^^^ 
while that for the non-degenerate case with p/f^^6is Up{s) — {^~Yi) ~ (y^-ra) iYi-Y2) 
with Yi = e-^l^ (-fiS + - ^) and Y^ = e'^l^ ("/^ - \l 7^ ~ ^) • inverse Laplace 
transform of these partial fractions with fractional power of s can be found only in the books 
of fractional calculus which give | (s-^-a) } ~ Sj=i (^^~^Et{j^—^, cfl) and | (^irz^j — 
El=i YL=i {tEt (U + m)v - 2, a<i) - [{j + m)v - 2] Et ((j + m)v - 1, a«)} with v = 

1, 1/2, 1/3, ... and q — 1, 2, Here the fractional exponential function Et{a, a) is de- 
fined as the fractional derivative of an ordinary exponential function Et{a,a) = J^e"* — 
r(a+i+i) ^i^h the derivative formula -^Et{a, a) = Et(a — /i, a). A linear combina- 
tion of the fractional exponential functions is a potential solution of the fractional differential 
equation because its functional form will not be changed after being performed derivative 
operator with fractional or integral order. By applying the inverse Laplace transforms to 
the partial fractions of Up{s), we can obtain the analytical solution of Eq. (9) as 

U,{t) = X [y,'E,{1/2, Y^) - Y^E,{l/2, Y^) + Y,e^^' - Y,e^^''] , (A5) 
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for the non-degenerate case with 7^ S; and 

(A6) 

for the degenerate case with (3/ = 6 which have been articulated in Eqs. (10) and (11) of 
the text. 
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FIG. 1: (Color online) (a) Dispersion relation near the band edge of a PBG material with directional 
dependent curvature is expressed by the effective-mass approximation with the edge frequency ujc- 
(b) Photon DOS D{uj) of the anisotropic PBG reservoir exhibiting cut-off photon mode below the 
edge frequency uJc- 
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FIG. 2: (Color online) Dynamics of single qubit in PBG reservoir, (a) Qubit's excited-state 
probability P{t) = \Up{t)\'^ with different detuning frequencies 5/(3 = (wio — ooc)/ P from the band 
edge frequency ojc- (b)von Neumann entropy S{t) = —Tr [p^{t)logp''{ty\ = —\j^logX+ — X-logX- 
of the non-Markovian {6/(3 < 0) and Markovian {5/13 = 2) systems. 
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(c) (d) 

FIG. 3: (Color online) Variation of concurrence dynamics with respect to the initial entangled 
degree a for the prepared states = a |01) + \/l — |io) and \^) = a |00) + Vl — (a) 
Concurrence for the prepared \ ^) states with qubit frequency lying inside the PBG (non-Markovian 
system, 5/(3 = —5) and (b) outside PBG (Markovian system, 6/f3 = 2); (c) Concurrence for the 
prepared |^') states with qubit frequency lying inside the PBG (non-Markovian system, 6/ j3 = —5) 
and (d) outside PBG (Markovian system, 6/(3 = 2) 
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FIG. 4: (Color online) Concurrence dynamics of the two-qubit system for the prepared Bell states 
1$)^ = (|01)±|10))/V2 in (a) and for 1^)^ = (|11) ± |00))/\/2 in (b). Steady values of concurrence 
in (a) and (b) with respect to the qubit detuning frequency 5/ (3 
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